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Abstract. A brief review of various methods to calculate radiative accelerations for stellar 
evolution and an analysis of their limitations are followed by applications to Pop I and Pop 
II stars. Recent applications to Horizontal Branch (HB) star evolution are also described. 
It is shown that models including atomic diffusion satisfy Schwarzschild's criterion on the 
interior side of the core boundary on the HB without the introduction of overshooting. Using 
stellar evolution models starting on the Main Sequence and calculated throughout evolution 
with atomic diffusion, radiative accelerations are shown to lead to abundance anomalies 
similar to those observed on the HB of M15. 
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1. Introduction 

In his book, Eddingto^ (1 19261) evaluated the 
equilibrium concentrations to be expected if 
atomic diffusion in the presence of differen- 
tial radiation pressure were efficient. He con- 
cluded that some heavy metals should then 
completely dominate the spectrum. Since light 
and heavy elements are present in stellar spec- 
tra, he concluded (§193, 196) that some mix- 
ing process made atomic diffusion inefficient. 
He suggested that it be meridional circula- 
tion (§199). His argument for the importance 
of meridional circulation was qualitative. In 
a later paper he evaluated quantitatively the 
meridional circulation velocity without further 
commenting on its efficiency in mixing stel- 
lar interiors dEddingtonl 1929). While it was 
realized in the 1940s that the presence of 
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the giant branch in clusters implied that stars 
could not be completely mixed, Eddington's 
argument seems to have had enough weight 
to prevent the proper calculation of the ef- 
fect of differential radiation pressure until 
the late 1960s even if Eddington had par- 
tially corrected his argumen t (see the 193 
correction on page xiii of lEddingtonl 119261) . 
Work however continued on gravitational set- 
tling in outer sola r layers dBiermannl Il937t 
IWasiutvnskil 119581: lAller & Chapmanl 1 1960h : 
the most impo rtant application was made to 
white dwarfs (ISchatzmanll 1 9451) . 

Instead of assuming that equilibrium 
concentrations were reached, one of us intro- 
duced the differential radiative acceleration 
ter m, g rar i, into the dif f usion velocity equation 
of lAller & Chapmanl (1 19601) and calculated 
anomalies to be e xpecte d in atmospheric 
regions (Mich audl Il970l) . Comparison to 
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observations of ApBp stars suggested that g ra( j 
plays a role in at least some stars. As large 
atomic data bases started becoming available 
in the 1980s, we realized that it would be 
possible to calculate g ra d throughout stellar 
interiors with good accuracy at the same time 
as the evolution proceeds. All composition 
changes can then be taken into account self 
consistently during evolution for the g m d and 
the Rosseland averaged opacity. It was first 
tried to make those calcul ations wit h TOPBase 
from the Opacity Pro j ect (lAlecian et al.ll 1993c 
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1995allbh but using those data we could not 



reproduce the concentration dependence of the 
Rosseland averaged OPAL op acities around 
the s olar center (see §5.1 of iTurcotte et al.1 
1 1998b and w e shifte d to using OPAL spectra 
dRicher et aT] [l998) to calculate g rad . The 
original spectra they had used to calculate the 
OPAL opacity tables dlglesias & Rogers! 1 1993. 
119961: iRogers & Iglesias! Il992l) were kindly 
made available to us by Iglesias and Rogers. 

In this review of g ra d in stellar evolution, we 
will first briefly describe how the calculations 
are carried out in stellar evolution codes (§[2]i 
and compare the values of g la d obtained using 
OP and OPAL data (§EQ- A few examples of 
the effect of g, a( j in Pop I and Pop II stars are 
described (§[3) and, finally, results obtained re- 
cently for HB stars are presented (§@]i showing 
that g la dS play a role in stellar evolution (§0. 

2. Calculations of radiative 
accelerations 



The particle transport equations are introduced 
int o a stan d ard ste l lar evolution code d e scibed 
in iProffittl (11994 iProffitt & Michaudl dl99ll) 
and VandenBerg (1985). For each spe cies one 
adds a force equation (eq. [18.1] of [Burgers 
119691) and a heat equation (eq. [18.2] of 
Burgers). Similar equations are written for 
electrons. It is generally assumed that each 
atomic species can be treated locally as being 
in an average state of ionization. One needs to 
know Zj, an appropriate mean of the number of 
lost electrons. 

The dominant term for each species con- 
tains g m d(A) - g as a factor, where g m &{A) is an 



appropriate average of the radiative accelera- 
tion over the states of ionization of element A. 
Over most of the stellar interior and for most 
of the evolution, it dominates transport even 
though the electric field, "diffusion" and ther- 
mal diffusion terms are also included in the cal- 
culations and are important for part of the evo- 
lution in some stars. 

Rosseland opacities, mean ionic charges, 
and mean radiative forces are calculated us- 
ing the same interpolation method, based on 
th e principle of corr e spond ing states in § 2.2 
of lRogers& Iglesias (1992); see in particular 
their equations (4) to (6). Interpolation weights 
are determined for a subset of the data grid, and 
used to interpolate locally all these variables. 

In first approximation, evaluating g rd a(A) 
amounts to calculating the fraction of the mo- 
mentum flux that each element absorbs from 
the photon flux. In stellar interiors, it takes the 
form: 



grad(A) = 



4nr 2 c X A 



r *„(A) 

Jo k u (total 



V(u)du (1) 



where most symbols have their usual meaning. 
The quantities K u {total) and k u (A) are respec- 
tively the total opacity and the contribution of 
element A to the total opacity at frequency u, 
with u and V(u) given by: 



hv 
kf 



and 



<P(u) = 



15 



4tt 4 ( e « - 1) 



2 ' 



(2) 



(3) 



The calculations of g la d(A) involve carrying 
out the integration over the 10 4 u values for 
each atomic species, A. This must be repeated 
at each mesh point (typically 2000 in our mod- 
els) and at each time step (typically 10 4 up to 
the HB). The Rosseland average opacity is also 
continuously recalculated making these calcu- 
lations fully self consistent with all composi- 
tion changes. This has some effect on stellar 
cores (see Fig. \9] and the last two paragraphs 
of §5.1 of ITurcotte et"ai]| 19981 for a discussion 
of the solar case). 
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Fig. 1. Correction factors used in evolution calculations except for Fe, see the text. Adapted from Fig. 6 of 
Rich er et all 19981 . 



2.1. Correction factors 

The g ra dS of equation ([TJ are corrected by two 
processes. First by taking into account the frac- 
tion of momentum that is given to the elec- 
tron in a photoioni zation, as first suggested 
in this context bv iMichaudl (l970) following 
Sommerfeldl (Il939t) . and second by calculating 
the effect of having many stages of ionization 
and not only an average one, which has come 
to be called the effect of redistribution. 

The sharing of the momentum between the 
ion and the ejected electron is caused by the 
process remembering the direction of the in- 
coming photon and emitting the electron with 
a distribution which is not exactly spherical. 
The correction to sphericity is small (~ v/c 
where v is the electron velocity) but for a 
given energy, the momentum of the electron 
being c/v times that of the photon, the cor- 
rection to momentum transfer is large. The 
effect is difficult to calculate except for the 
ground state of hydrogen but accurate evalu- 
ations now exist for a few other cases and con- 
firm t he generalization formulas which were 
used (Massacrier 1996; Massacrier & El-Murr 



ll996t lEl-Mudl999tlSeatonll995h . Using their 
results, iRicher et al.1 ( 19971) showed that for 



any shell n of an hydrogenic ion, one could use 
the simp le formula for t he effect of momentum 
sharing (lMichaudlll970l) : 



/io„(«) = 1 - K 1 - v)- 



(4) 



which has been used for all cases calculated up 
to now. 

When many states of ionization are 
present, one usually works, for convenience, 
with an average state of ionization. The sim- 
plest "averaging" is to use a single diffusion 
velocity for each atomic species and calculate 
it for the average Z of the atomic species, 



- = E,X(A,)Z,- 



(5) 



where the sum i is over the states of ioniza- 
tion of element A. The use of equation ([TJ is 
in line with this approach. However each term 
appearing in a diffusion equation has its own 
Z dependence and this approach is not always 
satisfactory. The g la( j term is the most sensi- 
tive to the averaging process since very often 
the largest g m d{Ai) is for a very unabundant ion 
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whose lines are desaturated. One could use an 
average of the form: 



< grad(A) >= 



£,D(A,)X(A,)g ra d(A,) 



(6) 



where g ra d(A,) would be due to the fraction of 
the photon flux momentum absorbed while the 
atomic species is in state i. However this as- 
sumes that one may consider each ionic state of 
the atomic species A as independent. In prac- 
tice this is a poor approximation since there 
are frequent ionizations and recombinations. 
Consequently one needs to compare the colli- 
sion and ionization times and take into account 
the various ionization routes. This rapidly be- 
comes complex and has been solve d in de- 
tail only for He dMichaud et alJll979h and Hg 
dProffittetal.111999^ with additional calcula- 
tions described in Gonzalez et al. (1995a0)for 
CNO. 

In the evolution calculations these two 
effects, momentum sharing between ion and 
electron, and averaging over the ions, are 
combined into a multiplicative factor (= 
[gmd (redistributed) / g ra d(not redistributed)] ), 
which we call a correction factor, applied to 
the result of equation ([]]). It was calculated 
for a number of species (see Fig. HI) using 
OP data for solar composition and tabulated 
as function of R e (= Ne/T 3 ) and \ogT. They 
are used for all compositions. Those shown 
in Figure Q] were calculated at \ogR e = 2.5 
which corresponds approximately to densities 
within main-sequence models. One notes 
that the correction factors are close to 1 .0 for 
5 < log T < 6. This is the T interval over 
which g lac |S play the main role in evolution 
calculations done up to now. 

For log T > 6 the main contribution to the 
correction factor comes from the sharing of 
momentum with the electron. For all species 
shown except Fe the correction close to the 
center comes from hydrogenic ions for which 
it is believed to be reasonnably accurate. It was 
decided not to apply t he correction t o Fe fol- 
lowing the argument of Seaton] dl997f) that f mn 
should not be applied to autoionization reso- 
nances which often dominate in non hydro- 
genic ions and especially for Fe. For log T < 5 
the correction can be large and should be con- 



sidered uncertain when they exceed a factor 
of 2. This usually occurs where the g m d(Ao) 
contributes a large fraction of the value ob- 
tained with equation ©. Since the diffusion 
coefficient of the neutral state is much larger 
than t hose of the ionized states (Mi chaud et alj 
1978) a small concentration of the neutral 
can lead to a large effect which however de- 
pends on the domin ant ionization processes 
dMichaud etaTIl 19791) . 



2.2. Radiative accelerations from OP 
data 

In so far as we know, all stellar evoluton cal- 
culations done with g la d have been done using 
OPAL data for equation ([1) and the correction 
factors described above calculated using OP 
data. However it is possible to do all calcula- 
tions with OP da ta using the spectra available 
from their server dSeatonll2005l : lMendoza et a"D 
2007). Those spectra use 10 4 frequency val- 
ues just as OPAL spectra but they are equally 
spaced in a modified frequency which takes 
th e local flux int ensity into account (see §2.2.1 
of ISeator3l2005l). I t is sim ilar to that proposed 
bv iLeBlanc et all (2000) and should have a 
similar effect on accuracy as discussed there. 
This mesh should lead to greater accuracy for 
log T < 5. However OP has two disadvantages. 
First it has fewer atomic species than OPAL 
and the second point is that their data base does 
not contain the Z which is needed to calculate 
the diffusion velocities. These would need to 
be calculated separately. 

A comparison of g ra( j calculated with OP 
atomic data to similar calculations with OPAL 
data is shown in Figure [2] The agreement 
is seen to be quite acceptable though one 
should remember the scale. There are differ- 
ences by factors of up to 2 which actually 
appear though they are rare. Similar compar- 



isons were made by Delahave & Pinso nneaultl 



d2005l) and Beaton! (|2007l) with data taken from 



iRicher et al.1 d 19981) with simi lar agreement . 
See for instance Figure 8 of ISeatonl d2007l) 
for Si. However those authors note greater 
disagreement with comparisons they m ake in 
solar models with iTurcotte et al.l (Il998l) . The 
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Fig. 2. Comparison of radiative accelerations cal- 
culated from OP (dotted line) and OPAL (solid l ine) 
atomic data. From Fig. [7] of Ri cher et al J 19981 . 



greater differences appear to come mainly 
from the correction factors discussed abovqj. 
These are always included in the stellar evo- 
lution calculations we made and so are in- 
cluded in the figures of Turcotte et al.1 ( 19981) 
where [Delahave & Pinsonn eaull (2005) and 
ISeatonl ( 2007 ) took data for their comparisons. 
If, for instan ce, one looks at Figure 1 1 of 
ISeatonl d2007l) . the difference between OP with 
mte and OPAL is mostly caused by correc- 
tion factors shown on Figure Q] above since 
"mte" contains only part of the corrections 
we include. The cor rection factor s are n ot in- 
cluded in Figure 1 oflRicher et al.1 d 19981) from 
whi ch figure jDe lahave & Pinsonneault d2005l) 
and ISeatonl d2007l) took the data for compari- 
son with that paper. 

The OP data has also been 
used to obta in semianalytic formu- 
las for grari (lAlecian & LeBlancl 120001: 
Alecian & LeBland 120021; lAlecian & Artru 
19901; iLeBlanc & Alecianl |2004j). These 
require much less computing power than 



1 There is how ever one error which we found 
in Figure 1 1 of iTurcotte et"aT] 1 19981 : the curve 
for C was plotted incorrectly. It does not cor- 
respond to the values used in the calculations. 
This error is responsible for the large st discrep- 
ancy found b y Delaha ve & Pinson neault ( 20051) an d 
Seator] d20071) . 



required by integrating over spectra. They 
have the further advantage of allowing an 
evaluation for species for which data is not 
available by using trends in spectroscopic 
properties. These are however less accurate 
than the detailed evaluation from OP or OPAL 
and do not allow to take into account the 
effect of individual concentration variations on 
Rosseland averaged opacity nor on g ra dS. 

3. Examples of the role of radiative 
accelerations in Pop I and II stars 

Stellar evolution calculations including the ef- 
fect of g ra dS have now been done for a large 
number of stars of both Pop I and Pop EL It is 
beyond the scope of this brief review to men- 
tion all effects of g ra dS. However we briefly de- 
scribe one effect found in Pop I stars and one 
in Pop II stars before giving some more details 
of recent results for HB stars (©. 

The largest structural effect of g ra dS in Pop 
I stars is the appearance of an Fe convection 
zone in all solar metallicit y stars more massiv e 
than about 1.5 M (see Ric hard et alj 1200 ll) . 
Iron contributes most to opacity at \ogT ~ 
5.3. Its radiative acceleration pushes iron from 
deeper in the envelope to the point where 
the g m d(Fe) starts decreasing. There Fe ac- 
cumulates during e volution (see Fig. [4] of 
iRichard et al.l 12001) approximately where it 
contributes most to Rosseland opacity. This 
leads to an important increase of the radia- 
tive gradient which, as the Fe abundance in- 
creases, becomes larger than the adiabatic gra- 
dient and an Fe convection zone appears. In 
a 1.5 M Q star this takes a significant fraction 
of the main-sequence life to occur. In stars of 
1.7 M G and more, this occurs very early in the 
main-sequence life. In stars with Z = 0.01, Fe 
convection zones start appearing at 1.3 M Q . 

The detailed treatment of the interaction of 
metals with H and He in diffusion processes 
and of the effect of concentration changes on 
Rosseland opacity has also shown that an ac- 
cumulation of metals occurs just outside con- 
vective cor es and causes semi- convection there 
(see §4 of IRichard et~ai]|200ll) 

The only requirement is for the region with 
logT ~ 5.3 to be stable enough for diffusion 
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Fig. 3. Difference between the radiative and adia- 
batic gradient around the He burning core boundary. 
In both panels, the lowest curve is at zero age HB 
and the uppermost curve is the final one calculated. 
In the calculations with diffusion (lower panel) the 
convective neutrality condition is automatically sat- 
isfied while it is not satisfied in the absence of diffu- 
sion. Overshooting is not included in these calcula- 
tions. 



processes to occur there. There is no adjustable 
parameter in those calculations. 

In Pop II stars, we were surprised to find 
that g ra dS and gravitational settling cause abun- 
dance anomalies by factors of 2-10 in a low 
metallicity cluster such as M92, much larger 
than those expected in higher metallicity clus- 
ters such as M5, M71 or 47 Tuc. Even in 
NGC 6397 which is only a factor of 2 more 
metal rich than M92, the effects are expected to 
be much smaller though they have appare ntly 
been seen in NGC 6397 dKorn et al.ll2006l) . In 
M92, ano malies might have been seen near the 
turnoff bv lKing et all (1 19981) but the noise was 
large enough to make this uncertain. This result 
would need confirmation. 



4. Diffusion in HB stars 

A 0.8 M model with Z = 10~ 4 was evolved 
from the ZAMS through the core expansion 
phase on the HB. Atomic diffusion was in- 
cluded throughout evolution and no adjustable 



parameter was involved. The results were com- 
pared to a similar model without diffusion 
dMichaud etap2007l) . 

It was found that atomic diffusion had little 
effect on age determinations using the HB lu- 
minosity but that there appeared two interest- 
ing effects on the structure of HB stars. After 
the discussion of instabilities at the junction of 
He burning cores and radiative zones presented 
by iPaczynskl dl970l) . it became accepted that 
the mixed C cores of HB stars were extended 
by overshooting or penetration to maintain 
convective neutrality at the boundary. It is cur- 
rently not possible to evaluate reasonably ac- 
curately the e fficiency of over shooting in HB 
stars. To quote lS weigarj (1 1 994|) , "Unfortunately 
the efficiency of convective overshooting un- 
der such circumstances is entirely unknown. 
Canonical HB theory assumes that the convec- 
tive overshooting is highly efficient...". But as 
may be seen from Figure [3] it was found that 
the presence of overshooting during the core 
expansion phase is made unnecessary by the 
presence of atomic diffusion. It was not shown 
that overshooting plays no role but rather that 
it is not necessary to invoke such a little under- 
stood process. 

Diffusion also causes an extension of the 
H burning shell into the He core. Through the 
outward diffusion of C from the core and the 
inward diffusion of H into the core there ap- 
pears an additional H burning region inside the 
He core which was calle d diffusion induced H 
burning (see Fig. [6] of iMichaud et al.ll2007l) . 
This leads to a slight increase in the ZAHB lu- 
minosity. 

4.1. Surface abundance anomalies 

Since HB stars are Pop II stars that just left 
the giant branch of globular clusters, they 
are all expected to have the same concentra- 
tion of_jn^tals 2 at least of those heavier than 
Al dGratton et alJl2004l) . The concentration of 
CNO and other relatively light species might 
show small va riations but Fe is not expected 
to be affected. IMichaud et al.l (Il983l) however 
suggested that g Ia &s should lead to overabun- 
dances of at least some metals in those stars 
where settling causes underabundances of He. 
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Fig. 4. Concentration of surface [Fe/H] exp ected 
in our HB models compared to observations of Behr 
( j2003|) for Ml 5. The continuous dark lines cover 
the interval from 5 to 30 Myr after zero age HB 
for a number of models whose mass is identified 
on the figure in M . The star marked with an ar- 
row has Vsin; ~ 16 km s" 1 . All other stars with 
Jeff > 11000K have V s ini <8km s" 1 . Adapted 
from lMichaud et alj|2008l . 

Glasp ev et al.l JT989) have confirmed the over- 
abundance of Fe in one star of one cluster but 
at the limit of detection and this observation re- 
quired confirmation. This predictio n has now 
been confirmed in many clusters dBehr et all 
1999tlMoehler et alj2000tlFabbian et alj2005t 
Pace et al.1 l2006h but in particular by iBehn 
(2003) for M15. Overabundances of Fe by fac- 
tors of 50-100 are seen in all HB stars with 
r e ff > 1 1500K while the cooler ones have the 
same Fe abundance as cluster's giants. 

In stellar evolution calculations, the sur- 
face concentrations depend on the exterior 
boundary conditions . In the calculations of 
Micha ud et al.l d2007l) . the simplest assumption 
was made, that of a mixed outer zone with- 
out any mass loss. The mixed mass was ad- 
justed to reproduce approximately the obser- 
vations of_Fe in one of the stars observed 
by iBehrl (120031) in Ml 5. The same model re- 
produced reasonably well the observations in 
other high T e ff stars of that cluster as may be 
se en in Figure ffl Further more, as may be seen 
in iMichaud et al.1 d2008l) the other anomalies 
are also reasonnably well reproduced. This is 
a striking confirmation of the role of g la dS in 
HB stars. 



5. Conclusions 

The availability of large atomic data bases has 
allowed, as described in §|2] to calculate stel- 
lar evolution models for Pop I and II stars 
up and past the giant branch including all ef- 
fects of atomic diffusion. They cause abun- 
dance anomalies not only in Ap star s, as orig- 
inally suggested by IMichaud! (11970 ). but also 
in HB stars of clusters Q4.lt and possibly in 
turnoff stars (^3|. They play an essential role in 
drivin g pulsations in sdB stars dFontaine et al.l 
120031) . the field analogue of HB stars. It is fur- 
thermore not only the surface region which is 
affected but 50% of the stellar radius and 10~ 3 
of its mass dRichard et alj2002t IMichaud et all 
120071). 

Eddingtonl d 1926b was right however in 
suggesting that competing processes also have 
a role to play. For instance, in Ml 5, rotation 
plays a role probably through meridional cir- 
culation. The arrow in Figure |4] shows the 
only star with T e s > 11500K which has no 
abundance anomaly. It is also the only rela- 
tively rapidly rotating star. lOuievv et al.l d2007l) 
have shown that meridional circulation ex- 
plained that the stars with T efi < 11500K 
have a normal abundance and not the 5x 
overabundance that Figure [4] would suggest. 
While atomic diffusion driven by g m ^s plays 
the main role in creating abundance anomalies 
on the HB, it has to compete with the effects 
of rotation just as in HgMn or AmFm stars 
dCharbonneau & Michaud! 199 lb . 
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